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ABSTRACT 

Finite element analysis of pipes buried in linear 
and non-linear elastic media has been carried out in the 
present work. 

Pipe-soil system is analysed as a soil-structure 
interaction problem. The pipe-interface-soil is taken as 
a single unit, where, the pipe is simulated by bar elements, 
interface by interface elements and soil by 4 — 1 3ST quadrilateral 
elements. Small deflection theory is used, lion-linear analysis 
has been carried out simulating stress-strain behaviour by 
hyperbolic law. To incorporate this in finite element analysis, 
incremental procedure has been adopted. 

The responses like springline thrust, crown pressure, 
crown deflection, springline deflection, maximum radial 
pressure and maximum tangential pressure have been studied 
for various values of parameters, such as radius-to-thickness 
of pipe ratio, modular ratio,. depth of oover-to-radius ratio 
and semi -load width-to -radius ratio for surface pressure. The 
responses have also been studied for the provision of soft 
trench material above the crown of pipe and non-homogene icy 
characterised by a linearly increasing modulus with depth for 
surface pressures. The responses of pipe have been studied 
for the variation of radius-to-thickness ratio, modular ratio 
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and depth of cover— to-radius ratio for gravity loads. The 
results obtained for surface pressure and gravity loads are 
compared with Burns elasticity solution. 



CHAPTER 1 


ETTRODTOTIOF 


1 , 1 General 

Buried pipes are used for v<rater supply as well as 
for drainage purposes under highways, embankments and 
railroads besides other applications. These are manufactured 
from a variety of materials in a variety of shapes and sizes, 
the most commonly used ones being corrugated metal and 
concrete pipe. They are subject to superimposed surface 
loads and loads due to the weight of the medium, in which 
they are buried. It has been observed that traditional 
methods used, generally overestimate the loads transmitted to 
the conduit. The methods assume the load distribution a priori 
on pipe through semi-empirical methods and hence do not 
properly represent the pipe-soil interaction. 

In a proper soil-structure interaction analysis it 
is recognised that the pipe and soil system act as a synergistic 
unit and the load distribution is determined in the course 
of solution. Development of powerful numerical techniques 
and high speed computers enable realistic modelling of 
interaction problems. The finite element method has proven 
to be a powerful tool in realistically modelling the soil- 
structure interaction behaviour and can account for non-linear 
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m aterial properties, complex geometries etc. 

1 .2 Problem Definition 

TTae object of this work is to develop and apply the 
finite element method to the analysis of interaction of a 
shallow buried circular pipe subjected to static superimposed 
loads and gravity loads. A planer-strain formulation of the 
interaction problem is attempted, simulating the pipe by 
bar elements, the medium by four noded quadrilateral elements 
and the interface behaviour by interface elements. Linear elastic 
and non>-linear elastic material models for soil behaviour 
are used and special cases of nonhomo gen© ity and effect of 
soft trench material on pipe responses are considered. 

1 . 3 Previous Studies 

Traditional empirical methods for design of flexible 
and rigid conduits are documented in Spangler(20) , Watkins(22). 

Most of the analytical and numerical methods 
available so far represent the soil system as a continuum and 
the pipe as a shell or a continuum. Solution techniques 
for the problem of soil -conduit system can be classified • 
either as analytical methods (using classical elasticity 
theory and shell theory to obtain ' exact' solucions) or 
numerical methods (using approximate techniques such as 
finite difference and finite element) . 
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Some of tiis earlier analytical solutions have been 
reported by Burns (5), Hoeg (12) and Luscher (18). luscher (18) 
assumed the soil surrounding the pipe to be a thick walled 
cylinder subjected to a uniform soil pressure around the 
soil cylinder. In his work theoretical expressions for the 
buckling resistance of the tube-soil system were developed 
based on the elastic theory and experimental work was carried 
out to provide the necessary correlation. 

Burns (5) and Hoeg (12) obtained closed form solution 
of the problem of medium-cylinder interaction analytic oily 
by developing static solution for an elastic thin cylinder 
embedded in a homogeneous, isotropic and elastic medium with 
surface loading. The no-slip condition and full-slip condition 
at cylinder-medium interface is considered. In the solutions 
for no-slip case, the radial and tangential displacements of 
the cylinder were equated to the corresponding displacements 
of the cylinder at the medium-cylinder interface. In the 
solutions for the full-slip case, the radial displacements of 
the cylinder and medium were equated, but the tangential shear 
stress was set to zero at the medium-cylinder interface. X 
uniform surface pressure 'p' was applied on two horizontal 
boundaries of the medium and a uniform pressure of ’Kp 1 was 
applied on the vertical boundaries of the medium (Fig. 1,1 ). 



Burns (5) and Ho eg (12) used about the ‘same approach to develop 
their solutionsHowever ,the resulting equations from Burns 
solution contain expressions for thrusts and moments in the 
pipe in addition to -expressions for displacements 

and stresses in the medium and hence are more useful to the 
designer. 

Several investigators, including Kay and Krizek (14) , 
Bar and Bates (7), and Lew (17) have examined Burns solution 
and proposed design methods based on this theory. 

Brown (3) used the finite element method to 
analyse the forces on rigid culverts under high fills. He 
used the incremental fill sequence technique for the analysis, 
but the soil medium was assumed to be elastic. Moreover, the 
interface boundary condition was treated only approximately 
by using triangular elements which were assigned very high 
stiffnesses compared to the fill. Brown et al (4) improved 
the representation of the culvert by using three overlapping 
triangular elements. Both rigid and flexible culverts were 
investigated to study the effects of slippage and no slippage. 

An and (2) studied the stress distribution around 
shallow buried rigid pipes using finite element method. The 
study was limited to elastic analysis and depth of embedment 
equal to one diameter of the pipe. The effect of width of 
imposed load, relative modulus ratio etc. on responses are 
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analysed and compared with closed form solutions. 

Abel et al (1) have studied the effect of pipe 
stiffness, depth of cover and effect of slippage at the 
interface for an elliptical pipe for surcharge loading. 

It was found that the effect of Poisson's ratio of pipe 
material and the soil is negligible and the slip at the 
interface is beneficial for establishing arching action, 
particularly for deeply buried pipes. 

Prakash et al (19) investigated the buried pipe 
under embankment using finite element method. The validity 
of assumptions in the Mar ston>~Sp angler Theory was studied. 

The parametric studies for height of fill-to -diameter ratio, 
diameter-to-thickness ratio and backfill properties were 
carried out. The effect of sequential construction by 
assuming linear and non-linear soil behaviour was also 
carried out. The effect of non-linear behaviour on responses 
was studied. 

Valliappan et al (21) analysed the buried pipe 
under surface pressure by assuming the medium to be elasto>- 
plastic, limited tension material. Conditions of slippage 
and no-slippage cases are studied by using interface elements. 
Two different yield criteria, namely, Von Mises and Drucker*s 
modified Von Mises, have been adopted for the elasto plastic 
idealisation of the medium. 
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Duncan (10) studied bhe “behaviour of long span 
mstal culverts by considering soil as a non-linear elastic 
material. He analysed metal arch culverts with and without 
curved stiffening ribs. Pull bond between soil and pipe was 
assumed. 

Katona et al ( 13 ) presented a unified computer 
methodology for the evaluation of buried culverts made of 
corrugated steel, aluminium and reinforced concrete and a 
class of plastic materials. The program uses Burns solution 
as level 1 and finite element method as level 2 and 3 of 
analyses. Analytical modelling features incremental construc- 
tion and non-linear behaviour. 

1.4 Objective and Scope of the Present Study 

In the present work, soil-structure interaction of 
shallow buried circular pipes subjected to static surface 
pressures and gravity loads has been investigated using 
finite element method. The variability and non-linearity of 
soil medium properties, complex geometries and interface 
conditions and simulation of construction sequence encountered 
in the solution of the above geotechnical problem warrant 
the use of a powerful and versatile numerical technique such 
as finite element method. 

The traditional methods of design of buried pipes 
are based on experience and tests and rely mainly on field 


f 
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experience and test results. They generally overestimate 
the design loads due to surface pressures and gravity loads. 
'Burns (5) analytical solution offers an accurate assessment 
of soil-structure interaction of circular "buried pipes but 
contains several simplifying assumptions. The medium is 
assumed to be elastic, homogeneous, isotropic and weightless 
and the pipe is assumed to be at infinite depth. As a result 
of the loading condition, the response of the cylinder is 
symmetrical about the vertical and horizontal diameters. 
Hence, the use of this theory to predict the response para- 
meters for limited burial depth is to be verified. Moreover, 
the surface static pressure is assumed to be acting upto 
infinity on horizontal planes. The gravity loads can again 
be taken only as an approximate equivalent overburden 
pressure equal to the weight of column of depth of burial. 

The interface behaviour can be either full -slip or no-slip. 

The finite clement method can model the interaction 
problem without any of the above limitations of Burns theory. 
Any realistic material model can he incorporated in the 
analysis. The weight of the medium, the limited depth of 
burial and the limited width of load or concentrated load 
oan he incorporated within the framework of the method. 

Any arbitrary shape of the pipe can also he studied. 

The parameters of 'interest in the study are spring- 
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line thrust, radial pressure, tangential pressure, ero'ra 
deflection and springLine deflection. The effect of ratios 

of modulii of pipe-lo-soil and radius— to— thickness of pipe on 
the above responses for different values of depth of bur i al- 
to -diameter (H/D = 1, H/D=2) ratios are studied, for static 
surface pressure of full width. The effect of the interface 
behaviour on response parameters are also studied. The 
response of pipe by finite element theory and Burns theory 
are compared. The material non>-linearity characterised by 
a simple non»-l inear curve is used. The interface nonlinearity 
characterised by a non-linear bond-slip relation is also 
used to model the interface behaviour realistically. 

The variations of response parameters due to the 
self weight of the materials for various ratios of modulii 
of pipe-to-soil and radius of p ip e-to -thickness are studied 
for two depths of "burial (H/D=1, H/D=2) « The approximation 
involved in the use of Burns solution for estimating the 
stresses due to self weight is brought out. 

The Boussinesq theory is used for estimation of 
vertical crown stress due to loads at surface in traditional 
methods. A finite element analysis where the inclusion is 
present can give more realistic value of vertical stress 
because the interaction effects can he considered. The 
vertical stresses are analysed for various load widths and 
are compared with the Boussinesq' s values. 
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fhe non-homo geneit y of soil characterised hy a 
linearly increasing modulus of elasticity with depth is 
analysed. It is a common construction practice to use a 
soft trench material above the pipe to reduce the load 
transferred to the pipe. Bevr relevant parameters like 
height of trenc li-to-rad ius of pipe, width of trentch-to- radius 
of pipe and modulus of trench mater ial -to ’-modulus of soil 
are studied for one surface pressure case to see the effect 
on response parameters. 

Chapter 2 describes the general principles used for 
the analysis in the present work. The closed form solution 
developed hy Bums are tabulated for the response parameters. 
The elements used for the discretization of the pipe, continuum 
and interface are described briefly. The technique used in 
simulating non-linear elastic material behaviour of medium 
and non-linear bond-slip characteristics of the interface 
are also described. 

In Chapter 3> the results obtained by Sums solution 
and finite element technique for the solution of pipe-soil 
interaction problem are discussed. The effect of modular ratio 
B /Be., radius -to-thic kness ratio (R/t) and effect of load width 
on the response parameters are given. The effect of material 
non-linearity on responses is also given. The variation of 
response parameters due to gravity loads for various modular 
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ratios and rad ius>-to --thickness ratios is analysed. The 
effect of non >-liomo gene i t y of medium, and presence of soft trench 
material on the relevant parameters are studied. 

Conclusions are summarized in Chapter 4 with suggestions 
for further research. 




CHAPTER 2 


METHODS OP SOLUTION 


2 . 1 General 

In this chapter, Sums closed form elasticity solution 

> 

and the finite element method of solution are briefly presented. 
Basic assumptions common to all solutions are: plane»-strain 
geometry and loading, small displacement theory. The elastic 
theory is more restrictive than the finite element method 
with regard to the scope of the boundaries, bonding between 
pipe and soil medium, non-linear material behaviour, non>~ 
homogeneities etc . 

2.2 Burns Solution 

This provides an exact solution for an elastic circular 
cylindrical shell embedded in an isotropic, homogeneous and 
infinite elastic medium with a surface pressure (Pig. 1.1). 
Shell theory is used for the conduit and continuum elastic 
theory is used for the surrounding infinite medium. The 
conduit>-medium interface is modelled only for two extreme 
boundary conditions: bonded interface, where both tangential 
and normal stresses are transmitted across the Interface 
and frictionless interface, where only normal forces are 
transmitted across the pipe. 



15 


'fable 2.1 gives tie elastic responses of the pipe 
such as radial and tangential soil pressure on pipe, radial 
and tangential displacements on pipe and thrust on the pipe. 
The responses are given for hot! bonded interface and friction*- 
less interface. 

2 , 3 Finite Element Solution 

A static formulation of finite element method can be 
equivalent! v derived from virtual work principle or from 
variational principle. This is well documented in references 
8 , 23. 

The virtual work in any structural system can be 
expressed as: 

f 5{2) VfdY = / 6 M T {Qj dS +• / 6 (uj^fjdV (2.1) 

Y S V 

where ifi - stress vector 

{oj = strain vector 

^uj = displacement vector 

= surface traction vector 

if? = body force vector 
< ->T 

ij = transpose of vector 
61 } = virtual vector 
S = surface of body 
Y = volume of body 

Strain energy term in virtual work statement (Equation (2.1)) 



H 


is written in terms of displacements by use of the consitutive 
relationships and strain-displac ement relationships; i, e. , 

1^3 = [C] 1 3 Y (2.2) 




Lb] W's 


(2.3) 


v/here 


fc] 

[b] 


= Constitutive matrix 
= Strain-displac ement matrix 


Using the above relationships, virtual energy can be written as 

/ cl7 = / [B? 6 iuj T TCI f3f W, dY (2.4) 

Y Y 

The actual form of 03 and Jb} is dependent on material and 
kinematic assumptions and will differ between pipe and soil 
models. Here, the finite element approximation is introduced 
by sub-dividing the domain V into a discrete self of elements 
interconnected at common nodal points such that continuity 
is maintained at all points on the boundary of the elements. 

The assemblage of elements and nodes is termed the finite 
element mesh. 


By inserting Equation (2.3) into Equation (2,4) and 
allowing the integration over the entire domain Y (to be 
represented by a summation of element integrations) , the 
virtual work equation yields the global equilibrium equation: 

fel \uj = [R] 

where [k! ~z \hj e 

Ki =BVj e 


(2.5) 
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Global stiffness matrix Jkj and load vector {rJ are tie 
summation of the element stiffness matrix and load vector 
respectively, given by: 

Me = l IS e dv e 

e 

. H e = •/* ®e H ds e + ■/* fB] T z f} dv e (2.7) 
s e e e 

The element stiffness matrix, J}c) , is the heart 
of the finite element formulation and provides versatility 
for modelling complex boundary value problems by assigning 
any group of elements, special material characteristics, 
loadings, and/or boundary condition. 

2.3.1 Discretisation 

Contin uum : The continuum with boundaries is given 

in Rig. 2.1. Because of symmetry one half has been taken 
for analysis in this study. The lateral boundary which is 
smooth rigid is fixed at a distance of 41 from the centre 
of pipe, beyond which the stresses are assumed to be free 
field. The bottom boundary which is rough rigid is fixed at 
one diameter below the bottom of the pipe (Reference 2) . 

The continuum is discretized by 4>-0ST quadrilateral 
elements (Rig. 2.1) , &. quadrilateral element consisting of 

four constant strain triangles with each node having two 
degrees of freedom is shown in Rig. 2.2 . Internal degrees 
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of freedom for 4’-C3T elements are eliminated using static 
condensation technique. The details of the stiffness 
matrix derivation are available in standard text books( 8 , 23 ) . 


Pip_e : The circular pipe is discretized by one dimen- 

sional axial force elements as shown in Pig. 2.3. The 


stiffness matrix is be given by 



( 2 . 8 ) 


This local stiffness matrix can be transformed to global 
stiffness matrix. 


Int erface : The interface of soil and pipe is 

simulated by one dimensional interface elements developed 
by Goodman, R.S. , et al (11). Interface elements allow 
consideration of two subassemblies meeting at a common 
interface (Pig. 2. 4 ), such that under loading the subassemblies 
may slip relative to each other with Coulomb friction; 
or separate or rebond. 

The application of these interface elements in the 

present problem is in the treatment of the pipe-soil interface. 

* 

The behaviour of the interface will be seldom linear, and 
so non-linear bond slip relation has been used. The theory 
is given in section 2.4.2. The general stiffness matrix 
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derivation can be given as follows. 


The interface element as shown in Pig. 2 . 4 , has 
length I, but zero width as the nodal pairs (1,4) and (2,3) 
initially have identical coordinates. It is shown in local 
coordinate system with the X axis along the length and 
origin at the centre. Each node has two degrees of freedom. 
The stored energy, 0 , in such an element is due to the 
applied force per unit length acting through the displacements 
and must be summed through the element 

-j +l/2 

0 - v f w, P. dx (2.9) 

•-L/2 1 1 

In matrix notation, 

<j +1/2 

0 = — / {wf T iPt dx (2.10) 

>-L/2 

where jwj = relative displacement vector 


can be expressed as 

r 

w 


{vi] 


top bottom } 
- w g ( 


w. 


top 


n 


w. 


bottom 


n 


l 

J 


( 2 . 11 ) 


and = the vector force per unit length is 

p?} = | Ps f (2.12) 


The vector \pjmay be expressedin terms of product of unit 
joint stiffness and displacement as 
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p =kw (2.13) 

>•« 

7 / he re lc = a diagonal matrix property expressing the jo in t 
stiffness per unit lengti in the normal and 


tangential direction and can be written as 


k = 


k c 


0 


O' 


k 

n 


Substituting Equation (2.13) in Equation (2.10) 

+V2 m 

/ - 1/2 f [Id {w} dx 

j-Xi/2 


(2.14) 


The displacement, {w},may be expressed in terns of nodal 
point displacements through a linear interpolation formula. 
If u^ and v^_ are the displacements in the tangential and 
normal directions respectively at any nodal point along 
the bottom of the joint element as shown in Pig. 2.^, then 

0 1+2x/li 0 

2x 


S bottom / 
J w_ 

< s 


.-bottom 

L w n 


1 

2 


i-l 

1 


0 


0 



and along the top of the joint element 

2x 


f top. 

topf’ 

w 

n 


l 


1 

T 


1+ 


o 


1+ 


o 

2 x 

“IT 


2x o 


0 


2x 

1 » 


-T7 



Thus Equation (2.11) become; 


w = 




-A 

0 


0 

-A 


-B 0 
0 


B 

0 


0 

B 


A 

0 


0 

A 
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in which A = 1 


2X 2 X 

-- - and B = 1+ — - 
L 1 


1 


Symbolic ally, ?wj = [hi |u[ 


(2.15) 


Substituting Equation (2.15) in Equation (2.14) yields 


1 +L / 2 m t _ 

</ = 5 - / 1/4 )u(. fDJ jkl fo) Au( dx 

-l/2 


(2.16) 


By evaluating the integral, 


2 he 


k = 


0 Ik 


2k n 0 


0 


-1 In. 


'1k n 0 


2 ho 0 -2 k, 


2 k. 


-2 k, 


-11% 0 


: hn 0 -21% 0 


0 1 k, 

21 % 0 
2 k 


0 

- 2 ]% 


•Iko o 


-Ik, 


-n 

0 

1S » 

0 

E.Pn 


(2.17) 


This local stiffness matrix is transformed to glohrj. stiffness 
matrix for further analysis. 

2 . 4 Non-Linear Analysis 


A constitutive law defines the relation between 
stress and strain on phenomenological observations. In 
order to represent this in finite element analysis the 
stress-strain matrix should be defined properly. 
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2.4.1 Soil Fon-linearicy 


Most soils seldom behave linearly, The widely used 
function for simulation of stress>-strain curves in finite 
element formulation is by Duncan et al (9) using Kondner 1 s 
(15,16) findings that the plot of stress vs strain in a 
triaxial compression test is very nearly a hyperbola. 

Pig. 2.5 illustrates sucb a relation, which can, be stated 
in equation form as 


or 


1 


CTj 


a +b-} S-j 


(2.18) 


or 


r ■ 

°1 

(crp^) 


a +^ 0 -] 


(2.19) 


It has been found by londner et al (15,16) that the values 
of * a 1 and 'hi' can be readily obtained if plot is made 
according to Equation (2,19) which is shown in Pig. 2.6. 
Since the compressive strength reaches its peak before the 
curve becomes asymptotic, failure stress difference ( cr -j > "^)f 
O’an be related to ^7 a factor R^; i.e. 


The tangent modulus at any level of stress can be expressed as 


®ta~ 


OCT 


= J?1 

(frft-aSj) 


= s^i - 


MdisL ) 2 


(2.21) 
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where E^ = Initial tangent modulus 

Equation 2.21 stows that E^is a function of stress level. 

Incremental procedure is applied to studv non-linear 
analysis. Load is applied in increments and at the begining 
of each new increment of loading an appropriate tangent 
modulus value is selected (Equation (2.21)) for each 
element on the basis of stress levels in the element. 


2.4.2 Interface Non-linearity 


Most analyses have been performed using one of the 
two limiting assumptions concerning the characteristj.cs of 
the soil-structure interface: (1) the interface is perfectly 
rough; or (2) the interface is perfectly smooth. But, the 
stress-strain behaviour of interface is seldom linear. 

This non-linearity of interface is simulated by assigning 
non-linear stress dependent properties of the soil and 
interface behaviour. The stress-strain curve is simulated 
by hyperbola as given by Clough et al (6). The shear 
stiffness value can be given as 


K 


st 





^tan 6 


) 


2 


( 2 . 22 ) 


where Kq = stiffness number 

6 = friction angle in degrees 

n = stiffness exponent 



22 


R.£ = failure ratio 

< 5 - = normal stress 

^ = interfacial shear stress and P, i =atmosyhe .‘ic pressure. 

ca> 

Kj,n,R^. are to be found from experimental direct shear test curves, 
ihe procedure to incorporate this non-linearity is same as ahe 
ona given in section 2.4.1. 

2.5 So il-Pip e-Interface Stiffhess Matrix 

The stiffness matrices of the sub -systems have 
been described. Per the complete system, the formation 
of the stiffness matrix will be the simple addition of 
the pipe stiffness, interface stiffness and soil stiffness. 

Using strain energy principle, 

Strain Energy of the = Strain energy of Strain energy of 

total system the pipe system the interface 

system 

+ Strain energy 

of the soil 

system 

So, the stiffness matrix of the total system can be expressed 

= fels + t K li + tBp (2 -U) 

= stiffness matrix of the total system 
= stiffness matrix of the soil sysrtem 
= stiffhess matrix of the interface system 
= stiffness matrix of the pipe system 


where [£j 

KL. 

[Bi 

M P 
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The final equilibrium equation can he written as 

[ k ] tuj = { r } ( 2.2 4 ) 

where, is the vector of unknown displacements 

-\R'- is the vector of nodal loads 

Substituting the boundary conditions in Equation (2.24) 
unknown displacements of the entire system are found out 
using Gauss elimination solution technique. Knowing the 
displacements of the entire system, forces in bar elements 
and stresses in the continuum are found out. The stresses 
in the continuum are calculated at the centroids of the 
quadr il at er al s . 

2.6 Soil Medium with linearly Varying Modulus of Hast ic it y , E 

Soil, particularly sands, have linearly increasing 
modulus with depth. To meet this requirements, a laterally 
homogeneous and depth dependent elastic medium is presented 
here. The variation in E is assumed as (figure 2.7) 

3(z) = S 0 (.1 + oc z) ( 2 . 25 ) 

in which 3 at z = 0 is S Q , a is a constant and the poisson's 
ratio , and the density of the medi'ioi are taken to be same 
throughout for computation pruposes. .Accordingly, the 
element stiffness matrices (Equation (2.6)) are given as 
follows, which is now a function of depth z. 

63 e = VeM 1 [°a!e Kl e 3 o (1+a z) ( 2 - 2 0 

0 
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in which, E is taken out of matrix jcj (now it is T0^| 
z is the depth of the centro id of the quadrilateral element 
from the surface, h is the thick.: ess of the quadrilateral 
element. If the new jkj is found out at the element level, 
assemblage of the elements can he carried out in the usual 
manner described earlier. 

2 . 7 Oomput er Pro grams 

The following four programshave been developed to 
analyse the plane-strain buried conduits subjected to 
plane strain static loading. In all the programs, if nodes 
fall on straight line and are equidistant, data for only 
the first and the last point of this group are needed, 
intermediate nodal point lata are automatically generated. 

2.7.1 LEPBP 

This is a linear Elastic Finite element program for 
the analysis of Buried Pipes. The applied finite element 
procedure is based on c ons tant>-strain-el ements . The use 
of triangular and quadrilateral elements is permitted for 
continuum, although the latter is used in this study. The 
program has been developed using the bar elements for the 
pipe. The program requires an userndef ined mesh. The 
program can take no slip condition between pipe and soil. 
The flow chart is given in Pig. 2.8. 
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Pro gram 
MAIN 

DATAIN 

ASMBL 

BARED 

QUAD 

GST 

STRESS 

BARSTR 

BAN SOL 


Details 

Galls subroutine DATAIN, ASMBL, BAIISOL and 
STRESS 

Reads number of nodes, number of elements, 
number of materials, numBer of surface tractions, 
material properties, nodal point data, Boundary 
condition code, element data, material type, 
surface tractions. 

Calls suBroutine QUAD and assemBles the 
quadrilateral stiffness matriooe and Hoads in 
tlie gloBal system, modifies the system equilibrium 
equation By introducing the Boundary conditions. 
Calculates the stiffness matrices of Bar elements 
and puts in the gloBal system. 

Divides quadrilateral into four triangular 
elements, calls suBroutine GST to compute the 
quadrilateral stiffness matrix. 

Calculates the triangular element stiffness matrix 
and adds the corresponding four matrices upto 
quadilateral stiffness matrix. 

Galls the suBroutine BARSTR, calculates the 

continuum stresses and strains, calculates the 

¥ 

radial pressures and tangential pressures. 

Calculates the axial displacements and axial 
forces in the pipe. 

Solves the modified system equilibrium equations 
for displacements By using Gauss elimination tedbniqtie. 
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2.7.2 TjSFSPI 

This is a Dinear ELastic Finite element program for the 
analysis of Buried Pipes with Interface simulation. The same 
elements, used in LEFHP, are used in the program for pipe and 
soil. Interface is simulated hy interface elements. This 
program requires an user>-def ined mesh. The program can take 
no-slip and full -slip condition between pipe and soil. 

Progra m Details 

MAH Galls subroutine DAIAIIT , ASMB1, BAN'SOL and 

STRESS’ 

DAT AIN Reads number of cards, number of elements, number 

of materials, number of surface tractions, .material 
properties, nodal point data, boundary condition 
code, element data, material type, surface tractions. 
ASMHEr Galls subroutine QUAD and assembles the quadrilateral 

stiffness matrices and loads In the global system, 
modifies the system equilibrium equations by 
introducing the boundary conditions. 

BARED Calculates the stiffness matrices of bar elements 

and puts them in the global system. 

JOINED Calculates the stiffness matrices of interface 

elements and puts them in the global system. 

Divides quadrilateral into four triangular elements, 
calls subroutine CST to compute the quadilateral 
stiffness matrix. 


QUAD 
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GST Calculates the triangular stiffness matrix and 

adds tie corresponding four matrices upto 
quadrilateral stiffness matrix. 

STRESS Galls the subroutines BARSTR and IH FACS, 

Calculates the continuum stresses and strains, 
calculates the radial pressures and tangential 
pressures. 

BARSTR Calculates the axial displacements and axial 

forces in the pipe. 

I1TFACS Calculates the relative displacements and shear 

and normal stresses. 

BANSOL Solves the modified system equilibrium equations 

for displacements by using Gauss elimination 
technique. 

2,7.3 FEFBPI 

This program is a Non-linear Elastic Finite element 
program for Buried Pipes with the Interface simulation. The 
non-linear material properties and non-linear interface 
behaviour are incorporated in this .program. The program 
requires an user-defined mesh. The program can take no-slip 
and full-slip condition. The elements used are same as in 
LEFBPI. The flow chart is shown in Fig. 2. 9 . 
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Program 

MAIU 

DATAHT 

I HIT I A 
ASMBL 

BARED 

JOINED 

BAITS01 

QUAD 

GST 


Detai ls 

Galls subroutine DAT AS, reads number of increments, 
calls subroutine IITITIA, solves for successive 
load increments by calling subroutines ASBffHL, 

BA. SOI and STRESS, controls incremental procedure. 

Reads number of nodes, number of elements, number 
of materials, number of surface tractions, material 
properties, nodal point data, boundary condition 
code, element data, material type, surface tractions. 
Initializes all the vectors. 

Calls subroutine QUAD and assembles the quadrilateral 
stiffness matrices and loads in the global system, 
modifies the system equilibrium equations by 
Introducing the boundary conditions. 

Calculates the stiffness matrices of bar elements 
and puts them in the global system. 

Calculates the stiffness matrices of interface 
elements and puts them in the global system. 

Solves the assembled banded stiffness matrix using 
Gauss elimination technique. 

Divides quadrilateral Into four triangular elements, 
calls subroutine GST to compute the quadrilateral 
stiffness matrix. 

Calculates the triangular element stiffness matrix and 
adds the corresponding four matrices upto quadrilateral 
stiffness matrix. 
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STRESS Galls subroutines CSTRES, BARSTR and INEACS, 

finds the total stresses for continuum. 

CSTRES Calculates the incremental strains and stresses 

for the continuum and total and incremental radial 
and tangential pressures. 

BARSTR. Calculates the incremental and total axial displare- 

ments and axial forces in the pipe. 

INEACS Calculates the incremental and total relative 

displacements and stresses. 

2.7.4 BURRS 

This program uses the BURES elasticity theory equations 
for the responses of the pipe. The program can take no~slip 
and full -si ip condition. 



Sable 2.1: E lasticity Soluti o ns for Pressure on Pipe and Pipe Responses 
for Bonded and Frictionless Interfaces (Burns(5)) 
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Fig. 2-8 Row chart for main 

(Program Lefbp and Lefbpi) 
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CHAPTER 3 


RESULTS AND DISCUSSIONS 


3 . 1 General 

The formulation described in the previous chapter 
is applied to solve the soil structure interaction of buried 
pipes. Linear analysis results' and non-linear analysis with the 
same initial tangent modulus of linear analysis results, are 
presented and compared with Burns solution. Ho -slip and full- 
slip cases have been analysed. Pull-slip condition in elastic 
analysis is achieved by using interface elements and giving 
zero shear stiffness value in interface elements. This is 
achieved in non-linear analysis by giving the following values 
in interface elements as given by Clough et al (6) on the 
basis of experimental results for sands ; 

6 = 0 . 1 ° 

Kj- « 1.0 

n = 0 .0 
R f = 1.0 

where 6 = friction angle in degrees 

K = stiffness number 

I 

n = stiffness exponent 
Rf = failure ratio 

These values used in non-linear analysis, give very low 



stiffness values. 


Ro-slip condition is achieved in elastic analysis 
without using interface elements or giving high values of 
shear stiffness for interface elements. This can he achieved 
in non>-linear analysis hy giving the following values suggested 
hy Clough et al (6) on the basis of experimental results for 
sands; 

6 = 35° 

K z = 35000 
n =1.0 
R f = 0.9 

These values used in non-linear analysis give very high 
shear stiffness values for interface elements. 

Separation between pipe and soil is prevented by 
giving a very high normal stiffness value in interface 

10 

elements. The normal stiffness value given is 0.7000 x 10 0f/l. 

Results are non-dimensionlised as follows: 

1 . The crown deflections and springline deflections 

are non>-diinensionalised by pxR/M g where p is the 
intensity of the surface pressure, M g is confined 
modulus equal to E g ( 1 /( 1+|) ( 1 -2| ) and R is the 

radius of pipe. 

The springline thrust is non-dimens ionalised by 


2 
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3. The crown pressure, maximum radial pressure and maximum 
tangential pressure are non-dimensionalised "by p. 

In the gravity load analysis, p will be equal toX^H 

where 

= unit weight of soil 
H = cover depth of pipe, 

3,2 Parametric Studies for Surface Pressure 

The effect of various parameters on the responses like 
springline thrust, crown pressure, crown displacement, springline 
displacement, maximum radial stress and maximum tangential stress 
are studied. The general variation of crown pressure, tangential 
pressure and springline thrust arc shown in Pigs. 3.1 to 3.3. 

The various parameters studied and responses obtained are 
given here. 

(a) Variation of responses with R/t for no-slip and 
full-slip condition for full load width. 

(b) Variation of responses with I^/B s for no-slip and 

full**slip condition for full load width. 

L 

The above studies are mainly concentrated to invest igs/te 
no -slip to full-slip factors (P) . The responses and Pec tors (p) , 
for different parameters involved are given in Tables 3,1 and 
3.2. Graphically responses are shown around the pipe for 
a typical case in Pigs. 3.4 and 3.5* 
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These studies reveal that no-slip to full-slip responses 
differ very little as can be seen from Pigs. 3.1 to 3.3 and 
Tables 3.1 and 3*2. So all other parametric studies are made 
for no-slip condition. 

3.2.1 Effect of Rad ias-to- Thickness Ratio 

The responses have been studied for variouhs radius-to- 
thickness (R/t) ratios for full load width. The interaction 
curves are investigated for the commonly used range of R/t=l5 
to R/t=40. Modular ratio taken for this study is 400 . The 
analyses have been done for H/D ratios 1 and 2. Table 3.3 
gives the details of the parametric study as shown at (A) . 

( i) Gro wn Deflection : The maximum diametrical deflection 
occurs at the pipe crown relative to invert. The variation of 
dimensionless crov/n deflection with R/t ratios has been plotted 
in Pig. 3. 6. The crown deflection increases with increasing R/t 
in linear elastic, non-linear elastic and Burns solution for 
both H/D ratios 1 and 2. But the large increase in the 
cylinder wall thickness changes the crown deflection only by. a 
small amount. The non-linear elastic deflections are around 
1.9 times of linear elastic deflections for the same initial 
tangent modulus. The linear elastic deflections do not 
differ with increasing H/D ratio for all R/t ratios. But 
the non-linear elastic deflections decrease with increasing 
H/D ratio for all R/t values. 



( ii) SpringL ine Deflection: The variation of dimensionless 
springline deflection with R/t is plotted in Fig. 3.7. It 
can he seen that springline deflection decreases with increasing 
R/t values. But the large increase in the cylinder wall 
thickness changes the springline deflection only hy a marginal 
amount. Burns solution agrees with linear elastic finite 
element solution except in lower values of R/t ratios. 

Non-linear analysis gives higher springline deflections than 
linear elastic analysis. The springLine deflection decreases 
with increasing H/D ratio in non-linear analysis hut the 
decrease is not noticeable in linear finite element analysis. 

(iii) Gro wn P r e s s ur e:, The dimensionless- crown pressure 
variation with R/t is plotted in Fig. 3.8. Grown pressure 
decreases with increasing R/t ratio. The decrease of cfown 
pressure in Burns solution is steep whereas it is gradual in 
linear and non-linear analysis. Non-linear analysis predicts 
higher values than linear elastic analysis. Crown pressure 
decreases with increasing H/B ratio in non-linear finite 
element analysis. This decrease is very small in linear 
finite element analysis. 

(iv) Spring Line Thrust; The variation of dimensionless 
thrust with R/t values is plotted in Fig. 3.9 (Thrust is 
maximum at the springline). Thrust decreases with increasing 
R/t ratio hut dimensionless values are more than 1 in all 
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the analyses. This indicates that the pipe draws load more 
than pxR. Thrust increases with increasing h/D ratio for 
all R/t values in non>-linear elastic analysis hut this is 
not seen in linear elastic analysis. Burns solution given 
higher value than finite element solutions. 

(v) Maximum Radial Pressur e: The maximum dimensionless radial 
pressure variation with R/t is shown in Rig. 3.10. The maximum 
radial pressure linearly decreases with increasing R/t. The 
large decrease in R/t ratio changes the maximum radial pressure 
only by a small amount. But in Burns solution it increases with 
increasing R/t ratio up to R/t=20 and heyond this ratio, it 

does not change. The maximum radial pressure slightly increases 
with increasing H/B ratio in non-linear analysis but its increase 
is unnoticeable in linear finite element analysis. 

(vi) Maximum T a ngent ia l_ Pressur e : The maximum dimensionless 
tangential pressure variation with R/t is shown in Rig. 3.11. 
Similar trends of maximum radial pressure are seen for 
tangential pressure also. 

3.2.2 Effec t of Modular Rat io 

The interaction effects have been studied for, various 
modular ratios (B c /S s ) . The study is carried out for K/t=20. 

The dimensionless curves have been drawn for linear elastic 
analysis, non-linear elastic analysis and Bums solution. H/B 
ratios 1 and 2 have been studied. Table 3.3 gives the details 
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of the various parameters involved (shown at B) . 

( i/e ♦ Ihe variation of dimensionless crown 

deflection with E c /B s is shown in Big. 3.12. The crown 
deflection decreases with increasing ^/S g ratio. Burns 
solution agrees well with linear elastic solution. But non- 
1 inear solution shows a deflection of 1,9 times of linear 
elastic deflection for H/D =1 for the same initial tangent 
modulus. The deflections do not change with increasing H/D 
ratio for all E c /E g ratio In linear elastic analysis. But 
the deflections decrease with increasing H/D ratio for all 
E c /E s ratio in non*-linear elastic analysis. 

( ii) Spri ngline Defle ction: The variation of dimensionless 

springline deflection with ^/E s ratio is shovm in Big. 3.(3. 

The springline deflection increases with increasing E /E ratio. 

c s 

It decreases with increasing H/D ratio in non-linear elastic 
analysis but this is not seen in linear elastic analysis 
Burns solution compares well with linear finite element 
analysis. Non-linear analysis predicts higher deflections. 

( i 11 ) pgpwn P ressure ; The crown pressure variation with 
E c /E s is shown in Big. 3.14. The crown pressure increases 
with increasing B^/E g ratio. But it remains constant beyond 
® 0 /E s ratio equal to 600. The crown pressure decreases with 
increasing H/D ratio for all E c / 3 g ratio in non-linear 
elastic analysis. But this decrease is not seen in linear 
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elastic analysis. Bums solution differs approximately 
by 0,1 p with linear elastic solution in predicting the 
crown pressure for all 3^/13 ratio , 

( i- v ) Springline Th ru st The dimensionless springline thrust 
variation with S c /E g is shown in Fig. 3.15. The springline 
thrust increases with increasing ^/S g ratio. The dimensionless 
thrust lies above 1,0 in all the analyses. Beyond 3_/S 0 equal 
to 600, the increase is nominal. The pipe draws more load 
with increasing TL/e_ ratio. Springline thrust also increases 
with increasing H/D ratio in non*-linear analysis but this 
increase is unnotic eable in linear analysis for all E^/s 
ratios. Burns solution predicts higher values than the 
finite element analysis. 

(v) Ma ximum Rad ial Pressure : The dimensionless maximum 
radial pressure variation with 3 i/33 is shown in Fig. 3.16. 

L* o 

The maximum radial pressure increases will increasing E c /E s 
ratio. It remains almost constant beyond S./E,, equal to 600. 
Maximum radial pressure increases with increasing H/D ratio 
in non-linear analysis for all E^/Eg ratio . But this change 
with H/D ratio in linear analysis is not noticeable. Burns 
solution gives higher value than finite element solution. 

( vi) M axim um Tange nt i al Press ure : The maximum tangential 
pressure variation with E^/E is shown in Fig. 3.17. The 
maximum tangential pressure does not change much with increasing 
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E c / E s ratio * It increases with increasing H/D ratio for 
all E c /E s ratio in non>-linear elastic analysis "but this 
increase is not seen in linear elastic analysis. Bur ns 
solution given higher value than finite element solution. 

3.2.3 Effect of Load Width 

The responses of pipe for various load widths have 
been studied. They are studied for a radius- to-thickness 
ratio of 20 and a modulus ratio of 400. H/D ratios 1 and 2 
have been studied. Both linear and non>-linear results are 
given. Table 3.3 gives the various semi ’-load width~to~ 
radius (b/R) ratios studied (shown at 0) . The variation of 
response around the pipe with change in semir-load width-to -radius 
ratio (b/il) is shown in Bigs. 3.18 to 3.20. 

( i) G rown Deflectio n: The dimensionless crown deflection 
variation with b/R is shown in Fig. 3.21. Crown deflection 
increases with b/R upto a value of 4 in linear and non-- 
linear elastic analysis and beyond this it starts decreasing 
for H/D =1 . The steepness of the curve decreases with 
increasing H/D ratio in both the analyses and beyond b/R=4 
the crown deflection tends to be constant, Fon-linear 
elastic analysis gives higher value than linear elastic 
analysis with the same initial tangent modulus. 

( ii) Springl ine Dsf 1 e o t ion : Similar trends of crown 
deflection are seen for springline deflection (Fig. 5.22) , 
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( iii) Crown Pressure : The dimensionless crown pressure 
variation with b/R is shown in Pig. 3.23. The crown pressure 
increases with b/R. HI the curves flatten beyond b/R=7. 

The crown pressure decreases with increasing H/D ratios for 
all b/R ratio . The elastic solution for vertical stress on 
the centre line for strip load on a semi- infinite mass is 
plotted. It is observed that for all b/R ratios considered 
the crown pressure obtained from finite element analysis is 
less than the one obtained by elastic solution for a strip 
load. The inclusion reduces the crown pressure by interaction. 

( iv) Sp ringline Th rust: The springline thrust variation 
with b/R is shown in Pig. 3.24. Similar trends as that of 
crown pressure have been seen in the variation of springline 
thrust also. Beyond b/R=4, the pipe starts drawing load 
more than pR in all the analyses for both H/D ratios 1 and 2. 

(v) Maximum Radial Pressure: The maximum radial pressure 
variation with b/R is shown in Pig. 3.25. The maximum radial 
pressure increases with b/R, There is not much of difference 
between linear elastic and non-linear elastic analyses in 
the computation of radial pressure distribution. The maximum 
radial pressure decreases with increasing H/D ratio for all 
h/R in both the analyse.s. 

(vi) Maximum Tangenti al P ressure ; The maximum tangential 
pressure increases with incx-easing b/R and beyond certain 



values of b/R, it starts decreasing (Rig. 3.26) « It readies 
peak at b/R=4 in non»linear analysis and at 3 in linear 
analysis for H/D =1 . 

3.2.4 Effect of Soft Trench Material: 

It is an usual practice to provide a soft trench 
material above the crovm as shown in Fig. 3.27. The purpose 
of this material is to develop positive arching. The 
influence of trench material width, height and material 
stiffness on the interaction effects have been studied. The 
responses are nondimensionalised using the responses of a 
homogeneous system, and so, response ratios less than 1.0 
indicate that pipe distress being relieved by the soft 
material. Table 3.4 gives the various parameters involved in 
this study. The response ratio can be defined as the 
response of system with soft trench material to the response 
of the system with homogeneous material. 

(i) Effect o f Tren ch Widt h: Thfe height of the trench 

material is taken 0.9 R. The study is made for trend' modulus 
equal to 1/10 E g and for E c /S g ratio equal to 400, H/D equal 
to 1 . 

The variations of springline thrust, maximum 
radial pressure and maximum tangential pressure with 
increasing trench. width are shown in Fig. 3.28. The springline 
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thrust and maximum radial pressure are reduced with increasing 
trench material width. They decrease by 12 percent at width 
1,1 R. Beyond width equal to 1.1R, the trench material 
does not help in reducing the stress response. But the 
maximum tangential pressure increases with increasing width. 

The tangential pressure distribution cannot be predicted 
because the distribution has got disturbed. 

The variations of crown pressure, crown deflection 
and springline deflection are shown in Fig. 3.29 , These 
responses decrease gradually upto a width R where the 
crown pressure, crown deflection and springline deflection 
get reduced by 25 percent, 15 percent and H percent 
respectively. Beyond this width, they show an increasing 
trend. The crown deflection will be more than the crowndef lection 
with homogeneous material for trench material width > 1.3R. 

The crown pressure and springline deflection approach a 
constant value beyond a width of 1.4R. 

( ii) Effect of T renc h Ma terial Hei ght: The responses axe 
studied for the variation of height of the trench material. 

The above study concludes that responses do not get reduced 
beyond width 1 .OR and so the width for the study is chosen 
to be equal to 1 R. The same radius >-tc -thickness ratio, 
modular ratio and K/D ratio have been studied. 
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The variations of springline thrust, maximum 
radial pressure and maximum tangential pressure 'are given in 
Pig* 3.30; The response ratios for springline thrust and 
maximum radial pressure drop down gradually to 0.88 at a 
height equal to 0.9R. Beyond this height the responses 
remain horizontal. Maximum tangential pressure does not 
get distressed with trench material. 

The variations of crown deflection, crown pressure 
and springline deflection with increasing height of trench 
material are shown in Pig. 3.31 • „ The crown pressure, 

crown deflection and springLine deflection drop down hy 
25 percent, 15 percent and 25 percent respectively at a height 
of 0.9R. iny further increase in the height of trench material 
does not give any improvement. 

( iii) Ef fect of Trench Material Mo d ulus The responses have 
been studied for various modulii of tronch material . Modulus 
of trench material is taken as a fraction of the modulus of the 
homogeneous soil. The same radius » to >-t hickness ratio and 
modulus ratio with reference the homogeneous soil have been 
taken. This study has been done for trench width equal to R 
and height equal to 0.9R. 

The variations of springline thrust, maximum radial 
pressure and maximum tangential pressure are given in 
Pig. 3.32, The response ratios for springline thrust and 
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maximum radial pressure drop down to 0.9 at 1/10* The 

response ratios increase with increase in the trench modulus 
"beyond this. The variation of tangential pressure does not 
show any definite trend. 

The response ratios of springline deflection, crown 
deflection and crown pressure reduce hy 25 percent, 18 percent 
and 25 percent fespectively at Ej-/E g equal to 1/10 (fig. 5. 33). 
The response ratios start increasing gradually to 1 with 
increasing trench modulus beyond this value. 

3.2.5 Uonhomogeneity Linearly Increasing Modulus of Elasticity 

of Soil with Depth 

Soils, particularly sands, have linearly increasing 
modulus with depth. 'The pipe responses are studied hy taking 
this variation into account. The radius-to>-thickness ratio, 
modular ratio and H/D ratio taken for this study are 20, 600 
and. 1 respectively. The responses are non-dimensionalised 
with reference to homogeneous soil>-pipe responses. 

The variation of modulus at any depth is given by 
Equation ( 2 . 2 $). Three slopes (oc) 0.05, 0.10 and 0.25 have been 
studied. Table 5*3 gives the parameters involved in this 
study, 

( i) Gr own L ef l^e ct ion a nd Sp ilngl in e Le flec tio n: Crown 
deflection and springline deflection variations with slope 
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are sho' n in Pigs. 3.34 and 3.35”. They decrease with 
increasing slope. The curves flatten at higher values of a } 
Crown deflection and springline deflection decrease with 
increasing H/D ratio for all slope . 

( ii) Crown Pre ssu re: As can he seen from Pig. 3.36 the 
crown pressures do not change much from the ratio 1 with 
increasing a for both H/D ratios 1 and 2. They are less 
than 1 at higher a values. 

( iii) Sprin gline Thrust, Maximu m Ra dial Pr essure and Ma ximum 

Tangen tial Pressure., It is apperent from Pig. 3.3 J 
that thrust decreases linearly with increasing oc values for 
both H/D ratios 1 and 2. it decreases with increasing H/D 
ratios for all values of oc. Similar trend has been found 
for maximum radial pressure variation(Fig,3.38). But the 
maximum tangential pressure decrease's linearly with cc but 
increases with increasing H/D ratios for all a values 
(Pig. 3.39) * 

3.3 Parametric Studies for Gravity loads 

The pipe responses due to gravity loads cannot 
be negLected, particularly at higher depths. In the present 
study, the responses due to gravity load have been studied for 
various modular ratios and radius ■‘-to-thickness ratios. 

As burns solution neglects the soil weight, an equivalent 
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load of X^E has heen applied as surface pressure for full 
width. Table 3.6 gives all the parameters involved in this 
study. H/D ratios 1 and 2 have been analysed. 

3.3.1 Effect of Jladius»-to>-Thiekness 

In this study, the modular ratio has been taken 
as 400. The responses are studied for various R/t ratio . 

( i) Grown Deflection : The dimensionless crown deflection 
variation with R/t is shown in Pig. 3.40. A value greater 
than 1.0 implies that the pipe soil system is softer in 
deflection resistance than a homogeneous soil field. The 
crown deflection increases with increasing P./t ratio for 
both H/D ratios 1 and 2. Burns solution gives slightly 

lower value than linear elastic analysis. The crown deflection 
decreases with increasing H/D ratios for all R/t ratios. 

(ii) Springli ne D eflec tion: The dimensionless springLine 
deflection variation with R/t is shown in Pig. 3.41. The 
springline deflection decreases with increasing R/t ratios. 

The springLine deflection also decreases with increasing H/D 
ratios for all R/t ratio . The change in crown deflection 
with larger variation in R/t values in Burns solution 
appears to be small. 

( iii) C rown Pressu re: The dimensionless crown pressure 

variation with increasing R/t ratios is shown in tag. 3. 42 

t,X' *■ n f 11 . . * ■ ■ 

nI k 62134' 
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The crown pressure linearly decreases with increasing 
R/t in linear elastic analysis. It also decreases with 
increasing H/D ratios for all R/t ratio . The crown 
pressure Decreases steeply with increasing R/t ratio in 
Burns solution. At R/t = 40 , the finite element solution 
gives 0.8 times of fot H/D=l whereas Burns solution gives 
only 0.65 times of 

( iv) Spri ngl in_ e_ T hru s t : The dimensionless springline 
thrust variation with R/t is shorn in Pig. 3.43. The 
variation of thrust is very flat and the dimensionless value 
is more than 1.0 for both H/D ratios 1 and 2. It reveals 
that the pipe drawn loads more than the soil column weight 
(negative arching) . Thrust decreases with increasing TT/D- 
ratio for all R/t ratio . The thrust value given by Bums 
solution is smaller than finite element solutions 

( v) Maximimi^Jladial^ Pr essure and Maximum Tajigential_ Pre_ssure_:^ 

Pigs. 3.44 and 3.45 show the variations of dimensionless 
maximum radial pressure and maximum tangential pressure with 
R/t. They decrease with increasing R/t ratios. They also 
decrease with increasing H/D ratios for all R/t ratio , In 
the case of Bums solution, the values remain constant 
beyond R/t=2G. 

3.3.2 Effect of Modular Ratio 

The responses of pipe have been studied for various 
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modular ratios. The R/t ratio used in this analysis is 20. 

( i) Grown D e_fl_ec ti on: The dimensionless variation of 

crown deflection with increasing 3 C /E S is plotted in Fig. 3. 46. 
Crown deflection decreases with increasing ratios for 

both H/D ratios 1 and 2. Burns solution gives slightly 
lower value than finite element solution. It decreases with 
increasing H/D ratios for all E c /S s ratio . 

(ii) Spring line Defl e ctio n: Fig. 3.47 gives the variation 

of dimensionless spring line deflection with ratios. 

The springLine deflection increases as F^/Eg ratio increases. 
Beyond I^/E s =400, the deflection does not change much with 
increasing %/s g ratio. The springLine deflection decreases 
with increasing H/B ratio for all \/~& 8 ratio . Bums solution 
differs with the finite element solution. 

(iii) Grown Press urej The cro?na pressure variation with 
E /E is shown in Fig. 3.48. The crown pressure increases 
gradually*' with increasing E c /S g ratio. Beyond B c /3 s =60O» 
the crown pressure remains constant. It decreases with 
increasing H/D ratio for all E 0 /E g ratio . Bums solution 
gives smaller value than finite element solution. 

(iv) Spr ingli ne Thrust: The springline thrust variation 
is shown in Fig. 3.. 49. It increases with increasing E c /E a 
ratio. Beyond %/3 s ratio equal to 400, it almost remains 
constant. The springLine deflection decreases with 
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increasing H/D ratio for all R /E_ ratio . Burns solution 

o s 

gives smaller values of springline thrust than finite element 
solution. 

( v) M aximum Radial Press u re an d Maxim um Tangential Pressure: 

The maximum radial pressure and tangential pressure 
variations are shown in Bigs. 5.50 and 5.51 . They increase 
with I^/E s Ratios, and beyond E c /S s ratio equal to 400, they 
remain almost the same. They decrease with increasing H/D 
ratio for all 3%/B s ratios. The finite element solution gives 
higher value than Burns solution in lower depths. 
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Table 3.1: J o Slip to Full Slip Resp onse Factors for the 
Variation o f R/t 

Jo slip to full slip factors (p) 


Parametric study R/t ratio for surface pressures 
(b/R= 8 - Pull load width) 

t = °* 3 . ^= 0.2 



1EF - linear Elastic Finite Element analysis 
KEEP- Ron-linear Elastic Finite Element Analysis 
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Table 3*2: 1 T 0 Slip to Pull Slip Response Factors For 

the Variation of S n /li„ 

ho slip to full slip factors (j?) 

Parametric study lS c /g g ratio for surface pressures 
(b/R =8 -*■ Full load width) 
i> s =0.3, V 0 - 2 




H/D 

b/R 

R/t 

Factors (F ) 







B c 

~300 

3 s 

E c 

- — =400 

v 

.. c=600 

ii, s 

a 

o 

-P 

Cl O 


1 

8 

20 
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S <D 

O i — 1 
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O CD 

& 

■ a 
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1 
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20 
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P © 

•H i — [ 
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1 

8 

20 
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1 
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"1 

CD 

H 1 

& 

1 

8 

20 
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H 

Pi 
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pres si 

Pq 

j| 

1 

8 

20 
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1 § 

H .. 

& 

$ 

1 

8 

20 
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0.9989 

1.003 0 

Spring 

thrust 

rp J 
fel 

1 

i 

; i 

i 

8 

20 
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IFF - linear Elastic Finite element Analysis 
NIFF - Ron-linear Elastic Finite Blemenc .Analysis 

















Table 3*3 J Parametric Studies for Surface Pre 
(b/R =8 -> Pull load width) 
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59 


Table 3.5: Study of Linearly Varying Modulus of So i l with 
Depth for Various Slope Indies 


iTon-homogeneous Soil-d! Cz) = S_ (l+otz) 


n /- 
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Table 3*6: Param etric Studies for Gravity. .L oads. 
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ig.3-29 Influence of soft trench material on responses for 
surface loading 
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ig * 3-31 Influence of soft trench 
surface loadina 
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Fig. 3 32 Influence of soft; t r one h motor i a I on r e sconces for 
' surface loading 1 
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1.3-34 Response ratio of crown deflection versus a 
for surface loading 



Fig 3-35 Response ratio of springline deflection versusa. 
for surface loading 
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.3-36 Response ratio of crown pressure versus <x 
for surface loading 
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Rg.3-37 Response ratio of springline thrust versus* 
for surface loading 
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Fig. 3-38 Response ratio of maximum radial pressure 
versus c* for surface loading 
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Kiy 3-39 Response ratio of maximum tangential 
pressure versus fo" sur face loading 
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Fig. 3-40 Dimensionless crown deflection versus 
radius -to- thickness ratio for 
gravity loading 
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Fiq.3-41 Dimensionless springline deflection versus 
radius -to- thickness ratio for gravity 
loading 
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Fig. 3-42 Dimensionless crown pressure versus 
radius -to- thickness ratio tor gravity 
loading 
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Fiq 3*43 Dimensionless springiine thrust versus 
radius -to- thick ness ratio for gravity 

loading 
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Fiq 3-44 Dimensionless maximum radial pressure 
versus radial- to - thickness ratio for 

gravity loading 
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.3-45 Dimensionless maximum tangential, pre: 
versus radius -to- thickness ratio tor 
' „■ gravity loading 
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Fig. 3-46 Dimensionless crown deflection vs modular 
ratio for gravity loading 
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Fig. 3-47 Dimensionless spring line deflection 

versus modular ratio for gravity loading 
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Fl 9‘ 3-4 8 Dimensionless crown pressure versus 
modular ratio for gravity loading 
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Fig 3 49 Dimensionless springline thrust versus 
modular ratio for gravity loading 
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Fig. 3-50 ■' "Oii^eft.sionlesis maxirmim radial pressure 
yer^ds-modutdr ratio for gravity loading 
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CHAPTER 4 


CONCLUSIONS AND RECOMMENDATIONS 

4.1 Conclusions 

The pipe-soil interaction problem bas been studied 
by using finite element method for a pipe buried in a 
medium. The effect of full -slip and no-slip has been 
studied by using interface elements. Linear analysis, and 
non-linear analysis with the same initial tangent modulus 
of linear analysis have been carried out. The non-linear 
interface behaviour has also been incorporated in the finite 
element analysis. Incremental procedure has been used to 
carry out the non-linear analysis of interface and material. 
The influences of the modular ratio of buried pipe and 
the surrounding medium, r ad iu s - to -thic hn e ss ratio of pipe 
and different surface load widths have been investigated. 

The pipe response for gravity loads has also been analysed. 

The effects of soft material above the crown and linearly 
varying modulus of elasticity of soil with depth for full 
width of surface pressure have been investigated. The 
following general conclusions can be drawn from these studies. 

1) Slip between pipe and surrounding soil does not 

seem to influence the responses. 



For full load width, Bures solution and linear 
elastic finite element solution are in good 
agreement in predicting the displacements. Crown 
deflections decrease with increasing modular ratio 
and increases with increasing radius -to -thickness 
ratio in all the analyses. ITon-1 inear analysis 
gives higher values of displacements. Radius-to- 
thickness and modulus ratios are not highly sensitive 
in predicting displacements. 

Crown pressure is less than the surface pressure 
for full load, width. But tho springline tnrust 
is more than the load above the- pipe > ie, pxR for 
full load width. R/t is more sensitive in predicting 
the crown pressure in Burns solution than the finite 
element solution within the practical range of this 
parameter . R/t and T^/Bg are not highly sensitive 
in predicting the springLine thrust. 

The design parameters such as thickness of the pipe , 
springline thrust etc. are not affected in linear 
analysis increasing R/T) ratios for full load width 
but the design parameters get affected in non- 
linear analysis. 

Crown deflection and spring line deflection 
increase with increasing load with. They start 
decreasing beyond b/R=4. The design parameters get 
affected in linear analysis' for increasing H/D 
ratio for all b/R < 8. 



Grown pressure increases with increasing b/R in 
"both linear and non>-linear analyses. The crown 
pressure gets affected with increasing H/D ratio 
for all b/R < 8 . The crown pressure predicted 
by finite element method is smaller t^an the 
traditional Boussinesq prediction (Rig. 3. 23). 

The crown pressure predicted by finite element 
method tends to coincide with Boussinesq' s value 
at higher b/R ratios. 

Spring line thrust increases with increasing b/R. 

The spring line thrust gets affected with increasing 
H/D ratio for b/R < 8 in both the analyses. Beyond 
b/R=4, the dimensionless value is greater than 1 in 
all the analyses ,*i.e. the pipe starts drawing load 
beyond b/R=4. 

In gravity load studies, all the responses decrease 
with increasing H/D ratio for all values of modular 
ratio and radius^-to-tbiclcncss ratio in linear 
elastic analysis which was not seen in the case v/hen 
only surface pressures are considered. 

In gravity loads study, Burns solution gives lower 
pressure value than finite element solution. 

The crown pressure is 20 percent less than the 
free field stress at R/t =20. 
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10) R/t is a sensitive parameter bat j3 c /ij g is an insensitive 
parameter m predicting responses due to gravity loads * 

1 1 ) Provision of soft material around the pipe improves the 
structural ciprity of the pipe. The optimum width, height 
and modulus of trench m iterial are R,0.9it and 1/10 of 
soil modulus respectively, 

12) Ifon-homogeneous analysis t iking linearly increasin; elastic 
modulus for the soil with depth gives lower response values 
than the ones predicted by homogeneous analysis, the 
decreases is linearly proportional to a, the slope angle . 
The responses decrease with increasin/ H/D ratio for all 
v-HLues of a. 

4.2 Re commenda Lions for Future rfork 

'the present work has shown the application of finite 
element method to study the responses of a pipe buried in soil 
by considering soil behaviour as linearily elastic -and non- 
lmearly elastic. Several topics for further research in this 
area are suggested below. 

the non-linear elastic property used in the work 
does not predict the f lilure of the soil. So, soil modelling 
can be improved by usin? better constitutive laws. 

The flexural stiffness of the pipe is not considered 
in this work. This is one of the important desi ;n parameters 
particularly for flexible culverts. This can be studied 
by simulating pipe by beam-column elements or shell elements. 

The constant-strain triangle elements have been 
used here. Higher order elements like 1ST, isoparametric 
elements may improve the accuracy of the responses. 



While the programs developed are qui be general for 
application to buried pipes of different shapes and sizes. 
Computations have only been carried out for buried circular 
pipes in this investigation. 'The studied can easily be 
extended to other commonly used pipe shapes. 
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